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ABSTRACT 


An orbiter and its payload substructure are linked together by five 
trunnion joints which have thirty degrees-of-freedom. Geometric com- 
patibility conditions require fourteen of the interface physical coor- 
dinates of the orbiter and payload be equal to each other and the 
remaining sixteen are free to have relative motions under Coulomb 
friction. This report presents the component modes synthesis method 
using fourteen inertia relief attachment modes for the formulation of 
the coupled system. The exact nonlinear friction function is derived 
based on the characteristics of the joints. Formulation is applicable 
to an orbiter that carries any number of payload substructures. 
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NOMENCLATURE 


A a* A — ’ a a 

V B — » b k 

V V D 

f f 

f 

s 

I 

k, k, K 


k a- V V k s 


V k p 


m, m, M 


N 


n 


0 


p k 

q 

s 


“b 


u j 

u 

s 

Greek Alphabets: 

A 

y, y 


inertia relief attachment modes, Eqs. (2.14), (2.21) 
free-free normal kept modes, Eqs. (2.12), (2.21) 
defined by Eq. (2.33) 
friction force 

applied force on substructure 
unit matrix 

stiffness, generalized stiffness, system stiffness 
matrices 

spring constants of joint shafts, Eq. (2.6) 

number of kept normal modes of orbiter and payload, 
respectively 

mass, generalized, system mass matrices, respectively 
normal force on nth free coordinates 
zero matrix 

generalized coordinate of attachment modes 
generalized coordinate of normal modes 
set of independent generalized coordinates, Eq. (2.20) 
transformation matrix of p to q 

coordinates of built-in end of joint shaft opposite 

to Uj 

free coordiantes of joint nodes 
constraint coordinates of joint nodes 
coordinates of interior nodes of substructure 

2 

diagonal matrix of m 

kinematic and static coefficients 
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e 


0 ) 

Superscripts: 
o, p 



angular displacement 

natural frequency of substructure 

orbiter and payload 
matrix transpose 
matrix inverse 
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I. INTRODUCTION 


1.1 Objective of This Research 

An orbiter and its payload substructure are linked together by five 
trunnion joints: two primary joints support load in the axial (x) and 

normal (z) directions, two stabilizing joints transfer forces in the 
normal direction and the keel joint carries load in the transverse (y) 
direction. The joint shafts of the payload and the support bearings of 
the orbiter are free to have motions relative to each other under Coulomb 
friction. The objective of this research is to make analytical formula- 
tion of this nonlinearly coupled system for dynamic response. 

1.2 Method of Approach 

Component modes synthesis has been accepted as the most efficient 
method for coupling substructures of large degrees-of-freedom for 
dynamic analysis in the last decade. The five nodes of the joints have 
thirty DOF's (3 translations and 3 rotations) of which the fourteen 
restrained coordinates are the interface coordinates and the remaining 
sixteen coordinates free to move under Coulomb friction are a subset of 
the interior coordinates. The orbiter and its payload are treated as 
two free-free substructures for which two sets of free-free normal modes 
are generated by their finite element models. These modes constitute 
the generalized coordinates. Fourteen inertia relief attachment modes 
for the interface coordinates form the dependent generalized coordinates. 

The friction forces acting on the free coordinates depend on whether 
the joints are in stuck state (no relative motion) or in a motion state. 
Any number of coordiantes may be stuck and the others in motion. Non- 
linear friction functions will be derived for the joints. A set of 
motion equations of the coupled system in terms of the generalized 
coordinates are obtained by using the geometric compatibility conditions 
and equilibrium of the joint interaction forces. 
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1.3 Component Modes Synthesis Method 

Since Hurty^ first proposed the method of coupling large sub- 
structures by component modes synthesis for dynamic analysis in 1965, 
numerous papers have been published, such as references 2, 3, and 4, 
for extension and modification of the original concept. Among the 
various versions of this method, the inertia relief attachment modes for 
unconstraint substructures and the residual attachment modes for con- 
straint substructures have proven to be the most accurate methods 
up-to-date. 

An outline and examples of the various approaches of the component 
modes synthesis can be found in chapter 19 of Craig's book^"^ . Because 
of the limited pages allowed for this report, the detail of this method 
will not be presented, but rather, the readers are referred to reference 
5. The formulation of the coordinates transformation matrix and the 
generalized mass and stiffness matrices in chapter 2 of this report are 
very similar to that given in Craig's book. Some of the Greek alphabets 
used in the equations in Craig's book are replaced by English alphabets 
for typing convenience. 

1.4 Description of Substructures 

Figure 1 shows the payload structure and its finite element model 
is shown in Figure 2. This structure is to be linked to an orbiter by 
five joints which can be seen from these figures. The joint shafts and 
the support bearings can be seen in Figure 3. Coordinates x, y, and z 
are axes along the body axis, transverse, and normal directions of the 
orbiter as shown in Figure 1. 

It is important for engineering designers to determine the loads on 
the payload structure and its interfacing structure accurately during 
the orbiter launching phase. The existing computing program has not 
been analytically formulated for dealing with frictional joints and is 
based on some simple assumptions on the joint friction forces. It is 
doubtful that such a simple model can produce accurate data for the 
designers. This is the motivation for this research. 
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(b) Keel joint. 

Figure 3. Trunnion joints of payload. 
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II. ANALYSIS 


2.1 Characteristics of Trunnion Joints 

An orbiter may carry one or more payload substructures and each of 
them is attached to the orbiter by five trunnion joints. As shown in 
Figure 3(a), the two primary joints support loads in x and z directions 
while the joint shafts are free to slide and rotate relative to the 
support bearings mounted on the orbiter. The two stabilizing joints 
shown in Figure 3(a) are designed to transfer load along z direction 
only and the joint shafts may have sliding and rotation motions in x 
and y directions relative to the orbiter. The keel joint which is 
located at the bottom of the payload structure as shown in Figure 3(b), 
has the same function as the stabilizing joint except the direction of 
y and z are interchanged. 

Consider that each interface node of the orbiter and payload has 

six degrees-of-freedom (3 translations and 3 rotations) which are denoted 

by u , u , u , 0 , 0 , and 0 . The geometric compatibility, or con- 
x y z x y z 

straint, conditions for the five joints are summarized in Table 1, An 

inertia relief attachment mode denoted by p is assigned for each con- 

a 

straint physical coordinate, u., as shown in Figure 4 and listed in the 
third column of Table 1 . The coordinates which have free motion rela- 
tive to the support, u^, are given in the last column. Discussion of 
inertia relief attachment modes will be given later. 

2.2 Coordinate Numbering for Joint Nodes 

Let u^ and u^ denote the physical coordinates of the interior and 
interface nodes, respectively. As shown in Table 1, among the 30 coor- 
dinates of the joint nodes there are 16 coordinates that have free 
motions (u^) which must be treated as a subset of u^ and the remaining 
14 coordinates form the set u.. These physical coordinates and the 
generalized coordinates for the attachment modes are arranged in 
sequence as shown in Table 2. 
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Table 1 . Characteristics of Trunnion Joints 


Joint 

u . Constraint 
J 

p Attachment Mode 
a 

u^ Free Motion 

Primary 

X C o 

II 

x - force 

u 

y 


0 p 

u = u r 
z z 

z - force 

0 

y 


e° = eP 

X X 

x - moment 



6° = eP 

z z 

z - moment 


Stabilizing 

o p 

u = u r 
z z 

z - force 

U 

X 


e° = e p 

z z 

z - moment 

u 

y 




0 

X 




0 

y 

Keel 

II 

c 

y - force 

U 

X 


e° = eP 
y y 

y - moment 

u 

z 




a 

V 

X 




0 

z 

Total Number 

14 

14 

16 


Superscripts "o" and "p" denote orbiter and payload respectively. 


Special notation is given to another set of physical coordinates 
which is a subset of u^, for the convenience of formulation of joint 
friction forces. These are the coordinates of the nodes of the joint 
shafts at the built-in end of the payload structure. The coordinates at 
the built-in end opposite to u^ at the free end is denoted by u^ as 
shown in Figure 4. 

In Figure 4 parts (a) and (b) depict the pairs of bearing and shaft 
of the two primary joints, parts (c) and (d) are for the stabilizing 
joints, and part (c) is for the keel joint. Nodes Bj through repre- 
sent nodes of the joint shaft at the built-in end. Other symbols are: 


XX IX- 7 





Table 2. Sequence of Joints Coordinates 


Sequence u^ (joint) 


1 

u (1) 

y 

2 

6 (1) 
y 

3 

u (2) 

y 

4 

9 (2) 

y 

5 

U x (3 > 

6 

u (3) 

y 

7 

9 (3) 

X 

8 

0 (3) 

y 

9 

U (4) 

X 

10 

u (4) 

y 

11 

9 (4) 

X 

12 

* 0 (4) 

y 

13 

% (=) 

14 

u 2 (5) 

15 

0 x < 5 > 

16 

V < 5) 


Uj (joint) p a (joint) 


U 

X 

CD 

f 

X 

(i) 

u 

z 

(i) 

f 

z 

(i) 

CD 

X 

(i) 

M 

X 

(i) 

0 

z 

(i) 

M 

z 

(i) 

u 

X 

(2) 

f 

X 

(2) 

u 

(2) 

f 

(2) 

z 

z 

CD 

X 

(2) 

M 

X 

(2) 

e 

(2) 

M 

(2) 

z 

z 


u 

z 

(3) 

f 

z 

(3) 

e 

z 

(3) 

M 

z 

(3) 

u 

z 

(4) 

f 

z 

(4) 

0 

z 

(4) 

M 

z 

(4) 

u 

(5) 

f 

(5) 

y 

y 


0 

y 

(5) 

M 

y 

(5) 


2d = distance between supports of the shaft in x direction 
L = length of joint shaft 
r =» radius of joint shaft 
r^ ■ radius of bearing support in x direction 

The above dimensions may have different values for different joints. 

The numbering system for all coordinates and forces are given according 
to Table 2. 
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2.3 Friction Force Acting On Free Coordinate 

2.3.1 Primary Trunnion Joint 1 

If the axial elastic force in the joint sh$f t is smaller than 
the axial friction resistance, there will be no sliding motion between 
the shaft and the bearing. With the aid of Figure 4(a), one may state 
the following. If 

Kf " u lbi k a < yN (2.1a) 

then, 

u lf " u if (stuck) (2.1b) 

f lf = “ f lf = k a (u lf “ u lb } (2.1c) 

where 

p = static coefficient of friction 
k = axial stiffness of shaft = 7 rr^E/L 
E = Young's modulus of elasticity 

N = pressure between shaft and bearing = [ (p ) 2 + (p ) 2 ]* 5 

i. a za 

Once the sliding starts, one may state that if 

U lf ^ U lf * (sliding) (2.2a) 

then 

f lf = “ f lf = k a (u lf ' u lb ) = yN (2.2b) 

When the elastic force in the joint shaft becomes smaller than the 
friction force, 

k al u lf “ u lbl < yN * (2.2c) 
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then sliding stops. 


Uj f * u,° (stuck) (2. 2d) 

where y is the kinematic coefficient of friction. 

In dealing with the coordinate ^f* two sets equations similar 
to Equations (2.1) and (2.2) may be written by replacing subscript "1" 
by "2", N by rN, and k by k = rrr^G/2L. The spring constant represents 
the shear stiffness of the shaft and G is the modulus of shear. The 
equations for the above two cases and the cases to be discussed for other 
joints can be represented by a simple figure as shown in Figure 5. Let 


x 


n 



u 


nb 




k x 
n n 


n 1,2,. ..,16 


(2.3) 


The function f ^ and the "stuck" and "sliding" conditions are shown in 
Figure 5. 


As x R increasing reaches point A where 


f c * k x = yr N 
nf n n n n 


(2.4) 


relative motion of the coordinate u^ starts; and then, the friction 
force drops immediately to point B where 


nf 


k x = yr N 
n n n n 


(2.5) 


If elastic force in the joint shaft decreasing passes point B, relative 
motion stops (stuck state). To the negative side points A' and B' are 
the equivalent points of A and B, respectively. It is important to 
note that the magnitude of N q is a function of the generalized coordinates 
P a which determine the magnitudes of A and B. Next, with the aid of 
Figure 4, the notations k, N, and r are derived for n equal to 1 through 
16. 
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2.3.2 General Function of Friction in Joints 

The following spring constants of the joint shaft are defined: 
2 

= irr /EIL (axial) 

k « irr^G/2L (shear) k = Trr^G/2L 

sx x (2 . 

k r * 4EI/L (end moment with zero deflection) 

3 

k^ ■ 12EI/L (end deflection with zero slope) 


One may define the following for Equations (2.3), (2.4), and (2.5): 


k n ‘ k s 

k n * k r 

k - k 
n sx 

k - k 
n rx 


n - 1,3,6,10,14 
n - 2,4,8,12,16 
n « 5,9,13 
n ■ 7,11,15 
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n = 1,3,6,10,14 


r 


n 


1 


r » r 
u 


n = 2,4,8,12,16 


r 

n 


r 


x 


n = 5,7,9,11,13,15 


N n - [(P U ) 2 + (P 2a ) 2 l n-1,2 

\ ' 'O-Sa ' 2 + " * 3>4 

»„ ’ [(p 9a )2 + <p 10a /d)2 ' p ' 5 - 7 

N „ * [(p lla )2 + (p 12a /d)2] ” * 9,11 


N = 
n 

[(p 13a )2 + 

2 15 

< p 14a /d > 1 

n = 13 

,15 



/ [<P 9a > 2 + (pHj) 2 ] 1 * 

if 


o 

U 5f 

N = 
n 

1 !<P 9a ) 2 + OW 2 ^ 

if 

U 5f 

o 

U 5f 


| [< p lla> 2 

+ (UNg ) 2 ]* 4 

if 

u 9 £ ** 

o 

U 9f 

N = 
n 

^< p lla> 2 

+ (Vg) 2 ]^ 

if 

U p = 

U 9f 

o 

U 9f 


(” p 13a> 2 

+ (pN ^) 2 ]* 5 

if 

u 13£ 

J ; O 

* u 13f 

N = 
n 

*t< p 13 a> 2 

+ ^b x 13^ ^ 

if 

u 13f 

o 

= u 13f 


n = 6,8 


n = 10,12 


n = 14,16 


2.4 Component Modes for Unrestraint Substructures 
2.4.1 Physical Coordinates 

Consider that the orbiter and its payload are two unre- 
straint substructures to be coupled together by the 14 interface coor- 
dinates described in Section 2.1. Now, let's begin the formulation from 
the equation of motion of a substructure, 

mu + ku = f (2.7) 
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The mass and stiffness matrices m and k are generated by finite element 
models, and f denotes external forces applied to the substructure. Note 
that the superscripts "o" for orbiter and ,, p 11 for payload will be used 
when the situation requires to identify the substructure, otherwise, 
expressions without these superscripts are applicable to both. The 
physical coordinates u are separated into two subsets: for interior 

nodes and for interface coordinates. As described in Section 2.2, 
the free coordinates of the joint nodes, u^, belongs to a subset of u^. 
Thus, u^ consists of u g , the physical coordinates of interior nodes, and 
u^. The m and k matrices must be rearranged and consistent with u in 
the form. 


u = [u u £ : u.] T = [u : u.] T 

o J. • • J x • J 

and so is the matrix f, 

f- it,t t : fji T - [e ± : f/ 

where 


( 2 . 8 ) 


(2.9) 


f g “ external load acting on interior nodes 
f £ = friction force acting on u £ 
f j = external load acting on = 0 

The dimension of these elements are as follows: 

u £ : 16x1 tu : 14x1 (for orbiter and payload) 

and dimension for u is arbitrary and different for orbiter and payload. 

s 

2.4.2 Free-Free Normal Modes 

Solving eigenvalues and eigenvectors of Equation (2.7) for 
free vibrations, one obtains six rigid-body modes (zero frequency) and 
free-free elastic normal modes as follows: 

(k - a)^m)b i = 0 (2.10) 
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with 


b^mb r “1 r = 1, 2, ...» 6 (2.11) 

b T mb = 1 
e e 

b T kb = co 2 . e = 1, 2 k (2.12) 

e e e 

where b and b are normalized nxl eigenvectors of rigid-body and elastic 
r e 

modes respectively, n is the dimension of u and k is the number of elastic 
inodes to be kept in the formulation (k <. n - 6). The modes higher than 
k are dropped to reduce the size of the coupled structural system. Now, 
denote the matrix of the eigenvectors of the kept elastic modes by 

Bj^ — (b ^ b2 ••• b^) (2.13) 

and let p^ be the generalized coordinate associated with the normal 

mode b^. Note that the orb iter and payload may have different values 

of k, i.e. , k ^ k . 

o p 

2.4.3 Inertia Relief Attachment Modes 

An attachment mode a is a set of physical coordiantes u due 

a 

to a unit force applied to one of the interface coordinates, u^, while 
forces on all the rest of the interface coordinates are zero. Let p 

a 

be a set of generalized coordinates associated with the attachment mode 

a ' (a * 1, 2 14), as shown in Figure 4. Let matrix denote the 

set of attachment modes, 

A a = ( a i a 2 ••• a i 4 ^ (2.14) 

Since the orbiter and its payload are unrestraint structures, the attach- 
ment modes are due to a self -equilibrated force system, 

f - f - mu 
e r 

which is the result of applied and rigid-body inertia forces. Hence, A a 
is called inertia relief attachment modes. For more detail the reader 
is referred to Reference 5, pages 479-487. 
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By definition, the inertia relief attachment modes may be obtained 
from the equation. 


A ■ G F 
a e a 


(2.15) 


where the elastic flexibility matrix may be formulated from 
G 


-1 T 
B A X B A 
e e ee e 


A ee " diag(l/m*) 


(2.16) 


and matrix F & represents unit force at each coordinate u^. 


F = 
a 


U ij 

L X jjJ 


0^ = ixj zero matrix 
Ijj = jxj unit matrix 


Thus, Equations (2.15), (2.16), and (2.17) result in 
A a ' B a£> 


(2.17) 


(2.18) 


Note that if only the first k elastic modes are kept for formulation of 
G g , and B g = as defined by Equation (2.13), and 


kk 


1/uk 


1/m. 


k 


(2.19) 


it can be seen that the contribution due to the neglected high frequency 
modes is small, if k is sufficiently large. 

It is easy to show that Equation (2.16) yields the flexibility. 

From Equation (2.12) 


A -1 - [B T k B ] _1 = B _1 k~ 1 (B T )~ 1 
ee e e e J e e v e' 


( 2 . 20 ) 


Premultiplying B e and postmult ip lying IT to both sides of the above 
equation leads to 


-1 T -1 
B A B - k - G 
e ee e e e 
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2.4.4 Coordinate Transformation 


It can be seen that the physical coordinate u is the result 


of the generalized cuuLuluafes and the 
interface forces p . Thus, 

Si 


.4.4 am « ■£ 4-V» a oann ynl ^ rrciA 

, W^VU \/JL. UtV 


u * (B k A a } 


V 

. P a. 


(2.21a) 


or in partitioned forms, 

irH'l 

L J jk ja a 


(2.21b) 


u 1 B . A 

s sk sa 

u r = B r1 A- 

f fk fa 

u . B .. A. 

J j Jk ja 


[::] 


(2.21c) 


2.5 Equation of Motion of Structural System 

2.5.1 Equation of Motion in Terms of Generalized Coordinates 

Substitution of the coordinate transformation and premulti- 
plication of the transpose of the transformation matrix given by Equation 
(2.21) to Equation (2.7) results in the equation of motion in terms of 
the generalized coordinates. 


“kk p k + kfck p k B sk f s + B fk f f + B jk P £ 


mp + k p = A ^f + + A.^p 

aa*a aa*a sa s fa f ja r a 


( 2 . 22 ) 


Making use of the orthogonality conditions given by Equations (2.11) and 
(2.12) and taking into account that A is a linear combination of the 
normal eigenvectors B^, it can be shown that the generalized mass and 
stiffness matrices in the above equation can be written in the following 
forms. 
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aa 


a a 


jk “kk kk D jk 


*kk 


kk 


^ ' E ak - 0 


E aa - A K - B jk ^Jk ■ A ja 


2.5.2 Equations of Constraints 


(2.24) 


Two sets of constraints on the interface coordiantes of the 
two substructures (orbiter and payload) that must be satisfied are 


Geometric Compatibility: 

v°. - u p - 0 
J J 

j = 1,2,.. .,14 

(2.25) 

Equilibrium of Interface: 

p° + p P = 0 
3 r 3 

a = 1 ,2, . . . , 14 

(2.26) 


Rewriting Equation (2.25) in terms of the generalized coordinates given 
by Equation (2.21) and putting Equations (2.25) and (2.26) into a single 
matrix equation, one obtains 


where 


Cp = 

0 



P = t p a p a 

p k p ki T 


C = 

0 cd 
<J^ 

1 

" A ja j B jk 

1 

VI 



1 i 0 

0 J 


[c 


1 

2 c 

aa { 


ak J 


(2.27) 


(2.28) 


(2.29) 


Let q denote the set of independent generalized coordinates p£ and p p 
and express the set of dependent generalized coordinates p° and p p in 

3 3 

terms of q. It gives 



-1 

r p k i 

„p 

= ’ C aa c ak « 

q - P 

L P a J 


L p k J 


where 


(2.30) 
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I - k,A * 


k = (A ° + A P ) ^ 
K 1 ^ ja ja ; 


(2.31) 


2.5.3 Transformation of p to q 

The complete set of the generalized coordinates p given by 
Equation (2.28) can be expressed in terms of the set of independent 
generalized coordinates q given by Equation (2.30) by the transformation. 


p = Sq 


(2.32) 


where matrix S can be readily obtained by using Equations (2.29) and 
(2.31) in the form. 


D = - k . B . j" 
o 1 jk 


d p ■ VjS 


D - [D 0 D p ] 


(2.33) 


2.6 Equation of Motion' of Coupled Structural System 

2.6.1 Equation of System in Terms of p 

Writing equations for orbiter and payload and putting 
together according to Equation (2.28), one obtains an equation of motion 
of the coupled structural system in the form. 


rap + kp = p s + p f + P a 


(2.34) 


where the elements are given in the following. The mass and stiffness 
matrices are 


(2.35) 
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The matrix P which denotes the external applied load, is 
s 

P , -‘C 0 B sk °1 T f s < 2 - 3 « 

and note that the external applied load on the payload is zero. The 
matrix P^ which represents friction load acting on the joint free coor- 
dinates u^, is 


[A. 


fa 


- A 


fa 


B 


fk 


-B Pl T f° 
"fk f f 


(2.37) 


and note that fP * — f^ • The term P a is due to the constraint forces 

acting on the interface coordiantes u.. With pP = - p°, one has 

J as. 


P = [A. 
a ja 


-A.P 

ja 


B ° 


-B P] T P° 

jk J a 


(2.38) 


2.6.2 Elimination of the Dependent Generalized Coordinates 

Substitution of Equation (2.32) into Equation (2.34) and pre- 

T 

multiplication of the resulting equation by S leads to equation of 
motion of the coupled structural system in terms of the independent 
generalized coordinates q in the form. 


M3 + Kq = Q s fg + Vf 

The elements in the above equation are given in the following, 
tern mass matrix may be written in the form. 


M = S T m S = 


where 


M M 
oo op 


, M M , 

L po ppj 


M oo ■ Hi + 

M PP - + Vjk 


where 


l l * < A ja + A ja> 


P'T 1 


M a M T = - (B °) T m.B P 
op po v jk 1 jk 

"l * k l ( "aa + S « )k l 


(2.39) 
The sys- 

(2.40) 
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The system stiffness matrix has the form, 
_ e'TjT o — 


IV — 


r*K 

K 1 


00 

op | 

(2.41) 

K 

K J 


po 

pp J 



where 


K oo • + (B 3k> T E l B Jk 


K = A. . 
pp kk 


0 + (.,!)' 


K = K 
op po 


- Vjk 


jk' 1 jk » 


k l = k l (k aa + k aa } k l 


The matrix Q , 
s 


r *o T + (b sJ )T ‘ 

<o T 


Q = SP = 
s s 


The matrix Q f , 


0- « S P. * 
'r r 


p sa 


D T (A f ° + A/) 1 + 
o v fa fa 


< B f£> T ] 

_ml 




(2.42) 


(2.43) 


The matrix product S L ? = 0. With the aid of Equations (2.29) through 

a 

(2.33), Equation (2.38) can be rewritten in the form, 

,T 


P = C 
a 


ra 


Then 


T 

SP 


a 


(CS) J 


[?] 


It can be seen from Equations (2.27) and (2.32) that 

CS » 0 and S T P = 0 
a 

2.6.3 Equation of Motion of Coupled Structural System 

Now, the formulation of equation of motion of the coupled 
orbiter and payload system. Equation (2.39), is completed. This equation 
is nonlinear as a result of the friction force being a nonlinear function 
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of the Independent variables. Therefore, the usual modal analysis of 
free vibrations can not be applied. However, one may examine three 
special cases from Equation (2.39). 

Case 1. Free vibration with frictionless joints (y ■ y ■ 0) 

Mq + Kq = 0 (2.44) 

Case 2. Forced response with frictionless joints 

Mq + Kq = Q Q f° (2.45) 

Case 3. Free vibration with frictional joints 

Mq + Kq = Q f f° (2.46) 

Both cases 1 and 2 are linear and can be treated by the common approach. 
The third case is nonlinear and its natural frequencies depend on initial 
conditions and amplitudes of motion. 

2.7 Orbiter Coupled with Two Payload Substructures 

Any number of payloads may be coupled with the orbiter at the same 
time. The modification of the formulation is simple and straightforward. 
To be specific, an orbiter with two payloads is illustrated in the 
following. 

2.7.1 Physical Coordinates Transformation 

The transformation given by Equation (2.21) is expanded as 

follows. 

(1) Coordinates of the orbiter 


s 


1 

td 

A 

sa^ 

A 1 

sa 2 





fl 


Bf i k 

A, 

£ i a i 

A- 

f l a 2 


' p k ‘ 


' B sk 

A 1 
sa 

rv 

f2 

= 

B f 2 k 

A. 

f 2 a l 

f 2 a 2 



= - 

B fk 

K 

fa 

L P a. 

jl 


V 

A. 

Jl a l 

A. 

Jl a 2 


- Pa 2- 


> 

A. 

J a _ 


j2. 


- B j 2 k 

A . 

j 2 a i 

A. 

J 2 a 2‘ 






(2.47) 
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where 


u ■ free coordinates of joints of payload 1 
It 


u. 


free coordinates of joints of payload 2 


u. * interface coordinates of joints of payload 1 


'1 


u, 


interface coordinates of joints of payload 2 


(2) Coordinates of payloads 


p f 1 

u 

s 


fiP s i 

0 

sa 

0 

P 2 

u 

s 


0 

sk 

0 

A P > 

sa 

p l 

u f 


B Pl 

B fk 

0 

fa 

0 

P 2 

u f 


0 

B ?2 

B fk 

0 


p l 
u . 

J 


£ 

0 

A P ‘ 

J 3 

0 

P 2 

u. 

_ J 


0 

p 2 

0 

A P2 

J a J 


■k (2.48) 

L 
a 

p 2 
3 

a 


where superscripts "p^' and "p 2 " denote payloads 1 and 2, respectively. 

2.7.2 The Generalized Corodinates 

The generalized coordinates defined by Equation (2.28) are 
increased and arranged in the form. 


r yJ u i 

P ■ b.j \ Pa Pa 


, o P 1 P 2,T 

q = [p v Pv Pv i 


P 2 ! 


’1 


2.T 


' F k F k 


P„ ] 


(2.49) 

(2.50) 


2.7.3 Other Symbols Used in the Formulation 

According to Equations (2.47) and (2.48), the following are 


defined: 
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0 


o A, 


also apply to A Jfl , A fa> B jk> B ja> and B fa< 


III. OUTLINE FOE NUMERICAL COMPUTATION 

3.1 Elements of Formulation 

Formulation of equation of motion of the orbiter and payload coupled 
structural system has been presented in Chapter II. Here, the essential 
elements for numerical computation will be discussed. Suppose that the 
finite element models for the orbiter and payload are available from 
NASTRAN, the first thing one must do is to identify the physical coor- 
dinates according to Table 2 and Equation (2.21). Then, do the following: 

(1) Make modal analysis of substructures to generate B (2.21) 

(2) Form inertia relief attachment modes A (2.18) 

(3) Compute the generalized mass and stiffness matrices: 

5 (2.23), k (2.24), M (2.40), K (2.41) 

(4) Formulate the loading functions: Q (2.42), Q f (2.43) 

s r 

(5) Program the friction force functions. Section 2.3. 

3.2 Friction Force in Terms of the Independent Coordinates 

The friction force function derived in Section 2.3 must be expressed 
in terms of the independent coordinate q for the solution of the equation 
of motion. The elastic deformation of the joint shaft given by Equa- 

tion (2.3) may be written by making use of Equations (2.21) and (2.33) 
in the form. 


p 

x - uK. 


-< 


(B fkPk + 


P D P 1 

'fa p a 


< B bk p k + 


» P _ R P 
B fk B bk 


<\a - A fa> D J 


(3.1) 


The normal force acting on the joints which causes friction Is a function 

of p which is given by Equation (2.33), 
a 

P° - Dq (3.2) 

a 
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As shown by Sections 2.2 and 2.3 and the above equations, the friction 
in the joints is a nonlinear function of q. 

3.3 Dimensions of Matrices 


In the following the dimensions of each symbol will be given in the 
parenthesis following the symbol. 


3.3.1 Physical Coordinates and Generalized Coordinates 

Let s^ and s denote the degrees-of -freedom of the inte 

op 

nodes, and k Q and k^ denote the number of kept normal nodes for the 
orbiter and payload, respectively. Then, 


u (s xl) 
s o 


u s (s p xl) 
Pk (k o xl) 

Pk (k p xl) 


q [< k o + k p) xl l 


u^ (16x1) 
uP (16x1) 


P a (14x1) 
pP (14x1) 


p [ (k Q + k p + 28)xl] 


u? (14x1) 


Uj (14x1) 


3.3.2 Transformation Matrix A and B 


Orbiter: 

i = s +16 
o 

j = 16, 

k = k , 
o 

f = 16, 

a = 14 

Payload: 

i = s +16 
P 

j = 16 

k = k , 
P 

f = 16, 

a = 14 


3.3.3 Other Matrices 


f s (s o xl) 


Q o [(k o + k p )xs o ] 


(16x1) 


Qf t(k Q + k )xl6] 
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VI. CONCLUSIONS 


1. The formulation is exact, i.e., there is no assumption or approxima- 
tion in dealing with the friction force in the joints. 

2. Formulation applies to structural systems with any number of pay- 
loads coupled with an orb iter. 

3. The equation of motion of the coupled orb iter and payload system can 
be used for the following cases: 

(a) Free vibration of a system with frictionless joints (linear 
problem) 

(b) Load response of a system with frictionless joints (linear) 

(c) Load response of a system with frictional joints (nonlinear) 

4. For free vibration of the system with frictional joints, the natural 
frequencies depend on initial conditions and amplitudes of vibration. 
The usual eigenvalue numerical procedure can not apply, since the 
system is nonlinear. 
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